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ABSTRACT 


We^prove  a  decay  estimate  for  the  steady  state  incompressible  Navier-Stokes  equations. 
The  estimate  describes  the  exponential  decay,  in  the  axial  direction  of  a  semi-infinite  tube, 
for  an  energy-type  functional  in  terms  of  the  perturbation  of  Poiseuille  flow,  provided 
that  the  Reynolds  number  does  not  exceed  a  critical  value,  for  which  we  exhibit  a  lower 
and  an  upper  bound.  Since  the  motion  is  considered  axi-symmetric  we  use  a  stream 
function  formulation,  and  the  results  are  similar  to  those  obtained  by  Horganv[5],xfor  a 
two-dimensional  channel  flow  problem.  For  the  Stokes  problem  our  estimate  for  the  rate 
of  decay  is  a  lower  bound  to  the  actual  rate  of  decay  which  is  obtained  from  an  asymptotic 
solution  to  the  Stokes  equations. 
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AN  EXPONENTIAL  DECAY  ESTIMATE  FOR  THE  STATIONARY 
PERTURBATION  OF  POISEUILLE  FLOW 
Gerardo  A.  Ache’ 


1.  Introduction 

In  this  paper  we  consider  the  steady  state  Navier-Stokes  equations  in  a  semi-infinite 
pipe  where  an  entrance  profile  is  specified  at  inflow  and  the  regularity  condition  that  the 
base  flow  becomes  Poiseuille  flow  is  specified  downstream  at  infinity.  Our  aim  is  to  estab¬ 
lish  a  spatial  decay  estimate,  in  the  axial  direction,  analogous  to  the  estimate  expressing 
Saint  Venant's  principle  in  elasticity  theory.  Early  results  regarding  the  mathematical  for¬ 
mulation  and  proof  of  decay  estimates  concerning  Saint  Venant's  principle  can  be  found  in 
the  works  of  Toupin  [14]  and  Knowles  [11]  for  linear  elasticity  theory  in  a  finite  cylinder. 
Later  this  type  of  result  was  extended  to  the  case  of  an  incompressible  viscous  flow  mo¬ 
tion,  by  Horgan  and  Wheeler  [8]  in  a  finite  cylinder  and  by  Horgan  [5]  for  two  dimensional 
incompressible  viscous  motion  in  a  semi-infinite  channel.  Horgan  and  Payne  [7]  have  es¬ 
tablished  a  decay  estimate  for  second-order  quasilinear  partial  differential  equations  in  a 
semi-infinite  strip. 

The  purpose  of  this  work  is  to  obtain  a  result  analogous  to  those  mentioned  above  for 
axi-symmetric  flow  in  a  semi-infinite  pipe,  i.e.,  we  prove  that  an  energy  functional  involving 
the  stationary  perturbation  from  Poiseuille  flow  decays  exponentially  fast  in  the 
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axial  direction,  provided  that  the  Reynolds  number  of  the  base  flow  does  not  exceed  a 
certain  quantity  Rc  (the  critical  Reynolds  number)  for  which  we  exhibit  a  lower  and  an 
upper  bound.  This  estimate  is  similar  to  the  one  presented  by  Horgan  in  [5]  and  it  is  proved 
using  a  stream  function  formulation  which  is  possible  since  the  fluid  motion  is  considered 
symmetric  about  the  axis.  In  [8],  Horgan  and  Wheek.  establish  a  decay  estimate  for  the 
three-dimensional  Navier-Stokes  equation  in  a  finite  pipe  of  arbitrary  cross  section  where 
Poiseuille  flow  was  specified  at  the  outflow  boundary.  In  both  results  !5]  and  [8],  there  is  a 
restriction  in  the  Reynolds  number  for  which  the  estimate  is  valid.  For  example,  in  [5  the 
estimate  is  valid  for  Reynolds  numbers  that  not  exceed  the  value  of  approximately  11.8, 
while  in  [8]  it  is  valid  for  Reynolds  numbers  less  than  approximately  4.8.  These  bounds 
for  the  Reynolds  numbers  actually  represent  a  lower  bound  for  Rc.  In  this  paper  we  have 
obtained  a  lower  bound  for  Rc  with  a  value  of  approximately  9.08,  also  we  have  obtained 
an  upper  bound  for  Rc  with  a  value  of  282.6.  This  last  number  arise  from  the  calculation 
of  the  supremum  of  the  quotient  of  two  integral  expressions  and  it  plays  an  important 
role  in  the  study  of  existence  of  solution  for  the  incompressible  Navier-Stokes  equations  in 
axi-symmetric  pipes  [3],  and  in  the  study  of  stability  of  Poiseuille  flow  [10],  [12]. 

Similar  to  :5j  and  [8j,  in  order  to  prove  the  decay  estimate  we  need  to  establish  some 
integral  inequalities,  for  which  we  provide  bounds  as  sharp  as  possible. 

In  section  6  we  compute  the  rate  of  decay  for  the  Stokes  problem,  for  the  stationary 
perturbation  of  Poiseuille  flow,  from  an  asymptotic  solution  to  the  Stokes  equations,  with 
a  value  of  approximately  4.47.  The  rate  of  decay  estimated  in  this  paper  has  a  value  of 


Wvvv 


Fv.*xr  *r  at  u"*  v*  v?  ir*  \r*  v*  v*' v^r*i 


approximately  1.79,  i.e.  the  rate  of  decay  estimated  in  this  paper  is  a  lower  bound  for  the 
actual  rate  of  decay. 


2.  The  Energy  Estimate 

Given  a  semi-infinite  tube  Rq  =  Co  x(0,  oo)  with  circular  cross  section  Co  of  radius  1.0, 
we  consider  the  steady  state  Navier-  Stokes  equations  in  cylindrical  coordinates  (r,  6,  2)  for 
a  viscous  incompressible  fluid  motion  which  is  symmetric  about  the  2-axis.  These  equations 
may  be  written  as 


duy  dux  dp  2  2 

u  i—  +  u2  — -  +  -r-  =  v(V  U!  -  ui/r*)  , 
or  dz  or 


du 2  du2  dp  o 

u  !—  +  u2—  -  —  =  u\-u2  , 
dr  dz  dz 


1  d  du2 


(2.1a) 


(2.16) 


(2.1c) 


where  u  =  (ui,0,U2)  is  the  velocity  field  and  p  is  the  pressure.  Here 


V2  =  1/r  d ! dr[rd j dr)  -d2 idz2  and  u  is  the  kinematic  viscosity. 

Associated  with  these  equations  we  consider  the  following  boundary  conditions.  At 


the  wall  of  the  tube 


«i(l,2)  -  «2(T  z)  =  0 


(2-1  d) 


and  at  the  center  u(0,2)  is  finite  for  all  2  >  0.  We  also  prescribe  an  entrance  profile  f  = 
(fi,f2)  such  that 


«i(r,0)  =  /i(r)  and  u2(r.O)  =  /2(r)  .  with  /i(l)  =  /2(l)  =  0.  (2.1c) 
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Finally  the  form  of  the  domain  leads  one  to  seek  a  solution  (u,  p)  of  (2.1)  such  that, 


u  — *  u°°  uniformly  in  r  as  z  —* >  oo  , 


(2-1/) 


where  u°°(r, z)  =  (0,0,  u(r))  is  the  Poiseuille  velocity  field,  i.e.  (u00^00)  is  a  solution  of 
the  steady  state  Navier-Stokes  equations  in  an  infinite  pipe  of  cross  section  Co  ,  and  for 
which 


I  udA  =  Q  >  0 

Jc0 


(2.2) 


Therefore  u  is  fully  described  by 


4(0  =  -  r*)  =  ^(1  -  '■*)  .  (2.3) 

41/  7T 

The  pressure  p00  is  such  that  its  gradient  is  given  by  Vp°°  =  (0,0.  -P)  ,  with  P  a  positive 
constant. 

We  define  a  perturbation  of  the  Poiseuille  flow  velocity  field  and  pressure  by 


vj  i  =  uj  ,  wi  -  u2  ~  u  ,  q  =  P  -  P°° 


(2.4) 


Then  the  pair  (w ,q)  =  ((w\,w?.),q)  satisfies  the  following  equations,  (perturbation  of 
Poiseuille  flow  equations) 


.  dw  i 

u~d7 


du>\ 

dr 


xv  2 


dwx  dq  2  tt’i 

=  ~  ~t) 


dr 

dw- 


.  du>2  dli  WU'2 

u~d7  *  d^U  l  *  Ul  "aT  ~  u’2~a 
i  a  au2 

;  (a?(-«',))  -  a7  =  ° 


dq 

dz 


=  uV2w- 


and  boundary  conditions 


U'i(l,-T)  =  u’2( l,c)  =  0  . 


(2.5a) 


(2.56) 


(2.5c) 


(2-5  d) 


Wl{r,0)  =  f\  (r)  ,  to2(r,0)  =  /2(r)  -  u(r) 


(2.5e) 


The  condition  (2.1/)  becomes 


w  — *  0  uniformly  in  r  as  z  — * ►  oo  . 


(2.5/) 


An  integration  of  (2.5c)  over  the  semi-infinite  tube  Rq  gives  us 


f 


2  7 r  j  w?dA  =  0 
J  Co 


(2.6) 


where  we  have  applied  the  divergence  theorem,  so  the  entry  profile  / 2(r)  must  satisfy 

2  t rf  f2{r)rdr  =  Q  (2.7) 

Jo 

We  now  proceed  to  describe  our  decay  estimate. 

Let  Rs  be  a  sub-region  of  R0  defined  by  =  C0  x  (s,  oo)  ,  similar  to  [5]  the  type  of 
decay  estimate  we  consider  involves  the  following  functional,  defined  for  0  <  s  <  oc. 


Vw  :  Vw dV  :  = 


dw 2  2  du'2  >2  u'r  .  ,  , 

i(a7>  '7*  ^  ' 


(2.8) 


Our  decay  estimate  is  stated  in  the  following  theorem,  which  will  be  proved  in  section  4. 
Theorem  2.1  Let  (u.p)  be  a  classical  solution  for  equations  (2.1a)  to  (2. If),  with  /2 
satisfying  (2.7).  Let  (w .q)  be  the  associated  perturbation  of  Poiseuille  flow  as  defined  in 
(2.5)  .  and  satisfying  the  following  conditions 


£(0)  <  oc  , 


(2.9a) 


f  Elc)dc  <  oc  for  all  s  0  . 


(2.96) 


w 2  is  bounded  , 


(2.9c) 


dw2  dwi  d2Wi 

— — ,  — — ,  -  -  are  bounded  uniformly  in  r  as  z  — *  oo 
oz  oz  oz2 


Then  there  exist  positive  constants  K ,  a,  v0  such  that 


E  (s)  <  KE(0)  exp(-as)  for  v  >  uq,  s  >  0. 


(2-9rf) 


(2.10) 


We  will  prove  this  theorem  using  a  stream  function  formulation.  The  proof  parallels 
the  result  of  Horgan  :o],  which  corresponds  to  the  case  of  two  dimensional  plane  flow.  The 
stream  function  formulation  is  discussed  in  the  next  section. 


3.  Stream  Function  Formulation 

The  main  difficulty  in  establishing  energy  decay  estimates  for  solutions  of  the  Navier- 
Stokes  equations  consists  of  eliminat  ing  the  perturbed  pressure  q  from  the  estimates.  In  8  . 
Horgan  and  Wheeler  found  a  device  useful  for  this  purpose.  They  consider  the  steady  state 
Navier-Stokes  equations  in  a  finite  pipe,  where  Poiseuille  flow  is  specified  as  an  outflow 
boundary  condition,  then  it  is  possible  to  eliminate  the  perturbed  pressure  by  the  use  of 
an  auxiliary  function.  However  their  approach  is  not  applicable  to  the  infinite  region  we 
are  considering  here. 

When  an  incompressible  flow  is  two-dimensional  or  axi-symmetric  one  can  eliminate 
the  pressure  by  introducing  a  stream  function.  It  is  then  possible  to  transform  the  per¬ 
turbation  of  Poiseuille  flow  equations  into  a  single  fourth-order  equation  from  which  the 
estimates  follows.  This  approach  has  been  considered  by  Horgan  5  .  for  the  case  of  two- 
dimensional  plane  flow.  Here  we  present  a  similar  result  for  the  axi-symmetric  case. 
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We  start  by  defining  the  following  stream  function  ip  > 


wi  =  - 


1  dip 
r  dz 


1  dip 
r  or 


(3.1) 


By  use  of  this  definition  the  incompressibility  condition  (2.5c)  is  automatically  sat¬ 
isfied.  Then  substituting  the  vector  field  w  by  (3.1)  in  equations  (2.5a)  and  (2.5b)  it  is 
possible  to  eliminate  the  perturbed  pressure  q  in  such  a  way  that  we  obtain  the  following 
fourth-order  equation, 


1  .dJ2iP  1  (dipdJ2xp 
-uW— +  — 


dip  dJ2ip  \  2  dtp  2 

dz  dr  /  r3  dz  ^ 


(3.2a) 


where 


ld_  ^  _  d_  \  d_  ^ 

r  dr  dz 2  dr  r  dr  dz 2 


We  may  integrate  ip  so  that 


(3.26) 


t P{r,z)=  /  w2(y,z)ydy  ,  (3.2c) 

Jo 

therefore,  by  integrating  (2.5c)  over  Rz  and  applying  the  divergence  theorem  we  have  that 
2tt  f*  w2(r,  z)rdr  =  0  ,  i.e.  ip  satisfies  the  following  boundary  conditions, 

ip  (0,z)  =  ip[l,z)  =  0  ,  --t/>(0,z)  =  -r^-(l,  z)  =  0  ,  (3. 2d) 

or  or 

~  (r, 0)  =  -r/i(r)  =  Fx(r)  ,  ip{r, 0)  =  [  (/2(y)  -  u{y))ydy  =  F2(r).  (3.2c) 

dz  Jo 

Condition  (2.5f)  becomes 


uniformlv  in  r  as  :  x. 


(3.2/) 


From  (3.2c)  and  assumption  (2.9c)  it  follows  that 


from  (2.9d)  we  have  that 


(3.2*7) 


1  (  d2rj)  d2ip  d3x p\  ,  j  j  .f  ,  . 

_  i  -  -  -  I  are  bounded  umiormly  m  r  as  z  — ►  oc  . 


H 


drdz  ’  dz 2  ’  dz 3  / 


(3.2h) 


Finally  the  energy  functional  (2.8)  in  terms  of  xp  is  , 


,  ,  0  f'  (,  d  .1  ,  d  1  dip  ^  .  d  .1 

d  1  dip  2  1  ,  . 


chb_ 

dr 


))2  - 


(3.3) 


We  now  proceed  to  establish  some  inequalities  which  are  useful  in  proving  Theorem 


2.1. 


Lemma  3.1 


Let  /  E  C![(C,  1)]  with  /( l)  =  0  ,  then 


f‘rdr  <  m  2  J  {~f)2rdr 


(3-4a) 


where  Hi  is  the  first  zero  of  the  Bessel  function  of  first  kind  and  order  zero  i.e.,  /r  i  =  2.4048. 


Proof 

The  function  /  can  be  expanded  as  f(r )  =  .4nJo(/rnr)  .  where  J0  denotes  the 

Bessel  function  of  first  kind  and  order  zero  and  the  /r„'s  satisfy  J o{nn)  =  0  ■  For  the 
functions  J0(ji„r)  we  have  the  following  two  orthogonality  relations. 

f  1  [  1  j  J 

J  Jo{nnr)Jo(vTnr)rdr  =  J  —  J0(finr)  —  J0(nmr)rdr  -  0  (m  n). 
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By  the  differential  equation  defining  Jo(r)  we  have  that 


_r  (r5r  J°(^"r))  =  ^  Jo{VnT)- 

Then  multiplying  this  last  expression  by  rj0(nnr)  and  integrating  by  parts  gives 

f  Jo  {Hnr)rdr  =  h\  f  {J'0{fj.„r)Y rdr  , 

Jo  Jo 

where  the  prime  on  Jo  indicates  the  derivative. 

Therefore  , 

/  f2rdr=YsAn  Jo(vnr)rdr  =  *22  ‘-2-  (  —  J0(f,nr))2rdr 

Jo  n_l  Jo  n=l  Mn  Jo  T 

-“:*L  • 

where  the  last  inequality  follows  since  n i  is  the  smallest  zero  of  J o(r). 

Lemma  3.2 

Let  /  G  C^O,  l)j  with  /( l)  =  0,  then  we  have 

/0  <s/)W)2  ’ 

where  o\  - 

Proof 

We  have  that 

/V)  =  -ij'mf/f  <  2 (J' /^y)'  '{J'(int)fdi)' 7  . 

so 

f4{r)r2dr  <  4  r2 f2{y)(^  f(t))2 dydtdr  . 


v  *.  ■,  v  *.  \  v  *.* •/  %  •;  •/ 


■v  ^  '  V  K  •"  *  * 


We  may  then  interchange  the  integration  in  t  and  r  ,  obtaining 


J  f4{r)r2dr  <  4  J  J  r2  f2{y)dydrdt 


Interchanging  the  integration  in  y  and  r  we  have  that 


r  1  r1  A  0  r 1  /-minty, t) 

J  }A{r)r2dr  <  4  J  (  —  /(f))  J  f2{y)  r2drdydt  = 

4  /  j  /2(y)§min(y’03^  <  ^  {^rf{r))2rdr  f2{‘Ardr. 


r2drdydt  = 


An  application  of  (3.4a)  gives  the  desired  result. 


Lemma  3.3 


Let  /' jr  £  Cx  '(0, 1) ;  with  /(l)  =  /'(l)  =  0  and  //r  — >  0  as  r  - *  0  ,  then 


/„  (Tr(lrTrn)2,ir  ' 


l  l~Tr{llf)'>2rdr 


(3.4c) 


(3.4<f) 


Proof 


An  integration  by  parts  gives, 


where  we  have  used  Schwarz's  inequality  and  (3.4a).  Therefore 


L  L  (7r('-r Jr n)~'Ti' 


For  the  second  part  of  the  lemma  we  apply  integration  by  parts,  Schwarz’s  inequality,  and 


(3.4a),  (3.46).  We  have 


('  r  .  ['  /3  <1  f  r  ,!  i  n  /  .  . 

L  =  /„  -'Tr>dr  =  j0  Mr 


1  f  2 


hence. 


Finally 


Then 


l  fMia'(j0  {T/;if^rdr)' 


Now  we  are  ready  to  prove  our  main  result. 


4.  Proof  of  Theorem  2.1 

To  prove  Theorem  2.1  we  use  an  approach  similar  to  that  used  by  Horgan  5  .  which 
consists  of  finding  an  integral  inequality  expression  which  leads  to  the  desired  result  (2.10). 
In  appendix  B  it  is  shown  that  E{s )  can  be  rewritten  as 


.  ,  f  f  xjj  a  f  d  rl.dy.,  1  .dv.n  1 

E{s)=2tt  /  -J4vdrdz-2n  —  -(— )2 - (-r-)‘ - 

J  Jr,  r  JCe  dzlr  dz  r  dr  r 


Cs 


-V 


d2vi 


dz2  J 


dr 


(4.1) 


and  since  y  satisfies  the  equation  (3.2a)  we  have  from  (4.1)  that 


E(s)  = 


271-  [  f  r  1  „  d  2  ,  1  dxP  dJ2v  dv  dJ2ii'] 

v  J  J  Rs  U  3:  r2  dr  dz  dz  dr  . 

„  f  d  \l.dxp  2  2  1  d2ri 

~2tt  - - (—  y  -  -  —  -  -V-r-^  dr 

JCa  dz  Lr  dz  r  dr  r  dz- 


\  ~  J2rl'\drdz 
r •*  dz  > 


Using  the  definition  of  J2  in  (3.2b)  and  ti  in  (2.4)  we  may  transform  some  of  the  terms 
in  (4.2)  as  follows, 


1 

-  iftlL  - - 

r  oz 


d  1.  2  d.udtpdv.  4 Q  dv  dv 

—  {-uvJ  v)  -  r-(--j-— ) - IT" 

dz  r  dr  r  dr  dz  tt  dr  dz 

1  d  (u  r  dv  j  ,  dxi'  .2)  1 
2^17 L(  1h’  ~  ( J  ’ 


ia,i  .dJ2v, 


(4.3a) 


1  (  dtp  dJ2tp  dip  dJ2tp  }  2  dtp  r2  1  5  ,  1 


^  <i6r  <9r 


K4- 


2  dr'  r2  dz 


1  d  1  ^dJ2v  d  1  2  t2  \ 

■ar)'fc(^r-7  '•'> 


(4.36) 


Then  substituting  these  two  expressions  (4.3a)  and  (4.36)  in  (4.2),  an  application  of 
the  Green’s  theorem,  the  boundary  conditions  (3. 2d)  and  the  regularity  conditions  (3.2/) 
and  (3.2 h)  transform  (4.2)  to 

E(s)  =  *“[(£)=  -  (frP]  -  ^v^r-v}dr 


12//  dJiaZidrdz)i  L  )J- i.^Ui 

tt  J  Jr,  dr  dz  J  JCs  dzlr'dz1  r  1  dr1  r 


t'lP2\dr\  '  ^4'4^ 


where  we  have  applied  integration  by  parts  on  the  term  .'  f  -[  i'2  j rJ2vdr. 


From  the  definition  of  E(s)  in  (2.3),  we  observe  that  the  derivative  of  E(s)  with  respect 
to  s  is  given  by, 

d  n  \  o  f  f(  d  l1  dvW2  i  d  ( 1  ■  /  d  i 1  d^\\2 

_£W  = -2, /{(-(--I)  -(s(-aj))  -  ar»  + 


,  d  ,1  ,  1  dV>2\  j 

(ai(r57l)  ~{7tY  \rdr 


1  d%\ 


Therefore  integrating  both  sides  of  (4.4)  between  the  limits  of  s  and  oc.  using  the 
expression  for  the  derivative  of  £(s)  in  (4.4)  we  may  write  an  identity  similar  to  one  given 


in  [5]  i.e., 


J  fOO 

±E(s)^4k2  E(‘)d‘  =  2  n(-h-I2~h)  , 

US  J  a 


where  , 


U  = 


7  [  hdtldli'))2  id(ldv))2  (d(ld^)Y 

l‘~ Jc,\ld~rYT;])  -{Jz{-t3;}>  _(s(;ar)) 

.  1  ^  2|  2  f  (1  dil'  2  2  1  d2Y' )  , 

+^ainrdr-*K  LUW  ~ - ;**?}*■  * 

4k2  f  f  <  u  ,  lur  Sti'.o  ,3v,2i  1  di>  t2  1  ,  . 

-t/  /J ;eJ  ^  jr  -  (57>  J  -  . 


(4.7a) 


(4.76) 


/3  =  - 


16Qk 


77/, 


dc  dv 
dr  dz 


drdzd£ 


(4.7c) 


and  where  k  is  a  constant  to  be  determined. 


Our  aim  is  to  give  an  estimate  of  the  integral  expressions  /].  /2,  and  /3,  by  use  of 
Schwarz’s  inequality  and  the  inequalities  stated  in  lemmas  2.2  to  2.4. 


1 

k 

§ 
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In  appendix  A  it  is  shown  that  for  k  <  k0,  where  k0  =  — )1//2  ,  the  integral 


term  l\  is  non-negative.  Also  in  appendix  A  it  is  shown  that  the  second  integral  expression 


/j  can  be  bounded  as 


I  2  <  4k ME{s)  , 


(4.8a) 


V  7T  7 r  '  -  1 


7 r  \  7T 


To  estimate  / 3  the  double  integral  term  in  (3.7 c)  can  be  bounded  using 


-H’i  • 


where  , 


It  =  sup 
<P 


f  f  do  do 

J  Jr,  dr  dzl 


drdz  )  /  |jd>;; 


(4.86) 


(4.10) 


(4.11) 


and  <{>  satisfying  the  boundary  conditions  (3.2 d).  Then  we  have  that  -)c  satisfies  the  fol¬ 


lowing  conditions. 


=  lo  >  0  , 


IS  < 


(4.12a) 


(4.126) 


It  >  1- 


(4.12c) 


where  -00  as  a  lower  subscript  stands  for  the  supremum  in  (4.10)  when  Rc  is  replaced  by 
an  infinite  cylinder  with  same  cross  section  as  . 


Condition  (4.12a)  is  obtained  by  a  shifting  argument  in  £.  (4.126)  follows  using  the 
arithmetic-geometric  inequality  ((A. 4)  appendix  A,  with  e  =  |)  and  (3.4a).  Finally, 
without  loss  of  generality  it  can  be  assumed  that  the  supremum  in  the  infinite  cylinder 
is  taken  over  functions  of  compact  support.  Therefore  we  can  approximate  7_oo 
arbitrarily  close  by  considering  only  functions  with  compact  support.  Then  for  any  <t>  of 
compact  support  in  the  infinite  cylinder  we  construct  a  function  v-  of  compact  support  in 
R0  by  a  translation  in  c  i.e.,  i'^(r,  z )  —  o (r,  z  —  k)  (for  some  k  >  0).  Therefore  the  quotient 
of  the  two  integral  terms  in  (4.10).  when  o  is  replaced  by  t? can  be  made  arbitrarily  close 
to  q.oo  ,  hence  70  >  7- 00  • 

Therefore  using  (4.12a)  we  obtain  the  following  bound  for  / 3  , 


8<?k27o 


/“sm  • 


(4.13) 


Using  these  bounds  for  /i,/2,/3,  the  identity  (4.6)  is  transformed  into  the  inequality 


r  00 

-E{s)  +  4k26  E{£)d$  <  4kME(s)  , 

J  S 


(4.14a) 


where 


6  =  1- 


(4.146) 


Notice  that  6  is  positive  for  values  of  the  viscosity  satisfying. 


4Q  7o 


(4.15) 


In  [5]  and  [8;,  Horgan  and  Horgan  and  Wheeler  have  shown  that  the  inequality 


rOC 

4 k2S  /  E(£)d‘  <  4 K.\1E(i 


(4.16) 


'j  'j*  */  O  '  •  V  .*  Jr 


leads  to  (2.10),  where  a  and  K  are  given  by 


a  =  2/c|\/ 


M2  +  6  - M 


and  K  = 


2yfW 


(y'M2  +  6  -  M) 


(4.17) 


respectively  and  M  and  S  given  by  (4.8b)  and  (4.14b)  respectively  (see,  e.g  [5]  pp.  371- 
372).  Therefore  by  taking  Vq  =  4Q~iq  n,  we  have  completed  the  proof  of  Theorem  2.1. 

Notice  that  using  (4.126)  we  can  obtain  the  following  upper  bound  for  Uq  i.e.,  i^o  < 
^s.  On  the  other  hand  using  (4.12c)  we  obtain  the  lower  bound.  u0  >  4 Q‘)-0o/Ti  • 

The  value  7_oc  was  computed  in  1907  by  Orr  (see  12  ,  pp.  157-159)  and  by  Amick 
(1978,  [4],  pp.  117-119).  They  computed  the  value  for  ~oc  of  approximately  1  359.82  . 
If  we  define  a  Reynolds  number  by  R  —  ^  ,  then  the  restriction  (4.15)  in  term  of  Reynolds 
number  becomes, 

R  <  Rc  =  ~  • 

47o 

Then,  on  using  the  upper  and  lower  bounds  for  7ch  Rc  lies  in  the  interval  9.084.282.6  . 

This  value  of  7_oo  plays  an  important  role  in  the  study  of  existence  of  solutions  for 
the  Navier-Stokes  equations  in  axi-symmetric  pipes  3  ,  and  in  the  study  of  stability  of 
Poiseuille  flow  [10],  !l2j. 


5.  Estimate  for  E(0) 

The  estimate  (2.10)  depends  on  the  value  £(0),  for  which  it  may  be  convenient  to 
provide  an  upper  bound.  As  pointed  out  in  5  and  8  .  the  estimation  for  E(0)  can 
be  given  in  term  of  a  functional  involving  an  appropriate  extension  of  the  solution  for 
the  Stokes  problem  in  Rq.  Such  an  estimation  will  introduce  another  restriction  to  the 


viscosity  or  Reynolds  number. 


The  Stokes  problem,  in  the  stream  function  formulation,  can  be  regarded  as  equation 
(3.2a)  for  which  the  viscosity  tends  to  infinity  i.e.,  if  rj  represents  the  stream  function  for 
the  Stokes  problem  then  we  have  that  77  satisfies, 


J  4r)  =  0  in  R0 , 


(5.1a) 


with  boundary  conditions 


V  ir,z)  =  ~i](r,z)  =  0  at  r  -  0. 1  , 


rj  (r,0)  =  F2{r)  ,  —  r?(r.O)  =  Fj(r)  , 


(5.16) 

(5.1c) 


and  the  regularity  conditions 


dt]  dt) 

( r?, —,—)—►  (0, 0. 0)  uniformly  in  r  as  z  cc 
dr  dz 


4  =  o(D 


(5. Id) 
(5.1c) 


and 


1  d2r)  d2r)  d3r) 

-  (  ■■  ■  ,  — r,  — — — )  are  bounded  uniformly  m  r  as  z  — >  oc 
r  ardz  dz *  dzJ 


(5.1/) 


We  define  an  energy  functional.  Eoc{s),  associated  with  r?  by  the  formula  (3.3)  with  v 
replaced  by  r)  where  the  infinity  as  lower  subscript  refers  for  the  infinite  viscosity.  From 
(4.8b),  (4.14b)  and  (4.17)  we  have  the  following  estimate  for  £cx>(5)  » 


E  <x{s)  <  2£oc(0)  exp(-2ics)  ,  s  >  0. 


(5.2) 


Also  it  is  possible  to  give  the  following  bound  for  E{ 0)  in  terms  of  ^^(O),  (the  derivation 
is  provided  in  appendix  B). 

E  (0)  <  .  (5  3) 

l  2  i/-  / 
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with  q  a  constant  given  by 


2^Qy  2Ex(0) 
nm  irn‘ 

•  In  order  that  (5.4)  be  well  defined  we  need  the  viscosity  be  sufficiently  large  so  that 

!  v  >  (2o) 1  2. 

t 

i 

» 

Using  a  modification  of  the  argument  in  13,  (see  appendix  C).  we  can  prove  that  for 
the  Stokes  problem  (5.1)  we  have  the  following  minimum  principle, 

£oo(0)  =  27rminSo(d>,d>)  ,  (5.5) 

<t> 

(the  bilinear  functional  S, ,  is  defined  in  (B.l)  appendix  B),  with  <?  satisfying  conditions 
(5.16)  to  (5.1e).  Since  for  any  function  6  satisfying  conditions  (5.16)  to  (5.1e)  we  have  that 
£oo(0)  <  2x Bq(4>i 0)  then  upper  bounds  for  £'oo(0),  and  consequently  for  £(0),  may  be 
obtained  on  constructing  a  suitable  auxiliary  function  6. 

6.  The  Rate  of  Decay  for  the  Stoke  Problem 

In  this  section  we  compute  the  rate  of  decay  associated  with  the  Stokes  equations. 
These  equations  are  obtained  by  setting  the  viscosity  equal  to  infinity  (or  the  Reynolds 
number  equal  to  zero)  in  equations  (2.1a)  to  (2.1c). 

Seeking  for  asymptotic  solutions  from  the  perturbation  of  Poiseuille  flow  equations 
(2.5a)-(2.5c)  in  the  form. 

w  i(r,z)  =  W'i(r)  exp(- Xz)  ,  (6.1a) 


(5.4) 


Q  =  ( 


w2(r,z)  =-  tt’2(r)exp(-  A;)  . 


(6.16) 


K  •' 


q(r,z)  =  <?(r)  exp(-Xz)  *  C  , 


(6.1c) 


with  C  an  arbitrary  constant,  yields  to  the  following  system  of  ordinary  differential  equa¬ 


tions 


1  =  AW, - -  , 


1  dll'. 


-  A  2 IV,  -  A  q  . 


dl  . 

dr  dr 


(6.2a) 


(6.26) 


(6.2c) 


with  boundary  conditions,  at  the  wall. 


w,(l)  =  W'2(l)  =0  . 


(6. 2d) 


In  order  to  compute  a  solution  to  this  system,  we  may  impose  the  following  symmetric 


condition  at  the  axis 


*V,(°)  =  ^(0)=0  . 


(6.2c) 


To  find  an  explicit  solution  to  this  system  differentiate  (6.2c)  with  respect  to  r,  mul¬ 
tiply  (6.2a)  by  A2  and  (6.26)  by  A.  then  add  the  resulting  expressions  and  using  (6.2c)  we 
obtain  the  following  differential  equation  in  q. 


d2i  .  \dl ..  xu  =  o 

dr -  r  dr 


This  equation  has  the  particular  solution. 


<?(r)  =  d,.(Ar) 


where  Jq  is  the  Bessel  function  of  first  kind  and  order  zero.  This  is  the  only  solution,  up  to 


a  multiplicative  constant,  which  is  finite  when  r  equal  to  zero.  Substituting  q{r )  in  (6.26) 
we  obtain  the  following  differential  equation  in  l-Vj, 

d2W2  1  dW2 


dr 2 


dr 


-r  A*H'2  =  —  AJc,(Ar)  , 


(6.5a) 


with  boundary  conditions, 


dW2 

dr 


(0)  =  ft'2(l)  =  0  . 


(6.56) 


This  two-point  boundary  value  problem  has  the  following  solution. 


W2(r)  =  bJ0(\r)~  -rJ^Xr)  , 


(6.6) 


wh°re  6  is  a  constant  to  be  determined  by  W2(l)  =  0  .  Using  (6.2a)  then  H']  is  given  as. 


^l(r)  =  ^^°(Ar)-(1-A6)Ji(Ar)  . 


(6.7) 


Since  the  solutions  (6.6)  and  (6.7)  need  to  satisfy  the  boundary  conditions  (6. 2d)  we  have 
that, 

A W i ( 1 )  =  -J0(A)  -  (1  -  A6)Ji(A)  =  0  .  (6.8a) 

2 

and 

W2(1)=U0{\)-  ^,(A)=0  .  (6.86) 

By  eliminating  6  from  these  two  equations  we  obtain, 

1  2  d0 ( A )  -  Jj,(A)  =  0  .  (6.9) 

2A  2 

To  solve  this  equation  we  observe  that  (6.9)  can  be  transformed,  after  some  algebraic 
manipulations,  into  the  following  non-linear  equation 

J  i(A)do(A)  -  ~(J?( A)  -  J2(A))  ---  0  .  (6.10) 

20 
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Then  the  stream  function  rj(r,z)  associated  to  the  equation  (5.1g),  has  an  asymptotic 
behavior  of  the  form 

r)(r,z)  ~  p(r)  exp(-aa)  ,  (6.11) 

where  a  is  the  real  part  of  A  satisfying  (6.10).  The  nonlinear  equation  (6.10)  can  be 
numerically  solved  using  Newton  method  1  .  the  rate  of  decay  a  is  given  by  the  real  part 
of  the  non-zero  solution  of  (6.10)  with  smallest  modulus.  We  compute  the  value  for  ct  of 
approximately  4.4663.  In  the  estimate  of  section  4  the  rate  of  decay  is  given  by  (4.17).  i.e.. 
for  the  Stokes  problem,  the  estimated  rate  of  decay  has  a  value  of  a  =  2/c o  approximately 
1.7906  which  differs  from  the  computed  rate  from  (6.9)  by  a  factor  of  approximately  2.5. 
That  is,  the  estimated  rate  of  decay  given  by  (4.17)  is  a  lower  bound  foi  the  actual  rate  of 
decay. 

7.  Conclusion 

An  energy  estimate  of  the  form  (2.10)  has  been  presented  for  the  Navier-Stokes  equa¬ 
tions  in  a  semi-infinite  pipe  for  a  motion  which  is  axi-symmetric.  The  main  result  is  stated 
in  Theorem  2.1.  The  estimate  depends  on  the  values  .£(0)  and  /c,  where  £(0)  is  assumed 
to  be  finite,  this  assumption  is  typical  in  the  study  of  existence  for  solutions  of  the  Navier- 
Stokes  equations  3  .  The  value  of  /c0  is  given  by  k  approximately  .8956.  For  the  Stokes 
problem  k  represents  a  lower  bound  for  the  actual  decay  rate  w’hich  is  given  by  the  real 
part  of  the  solution,  with  smallest  modulus,  of  the  nonlinear  equation  (6.10).  Finally  for 
this  proof  of  the  exponential  decay  estimate  it  is  required  that  the  Reynolds  number  does 
not  exceed  a  critical  value  Rc  which  lies  in  the  interval  9.08.282.6  . 
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APPENDIX  A 


In  this  appendix  we  will  derive  bounds  for  7j,  72  which  are  define  by  (-4.7a)  and  (4.7 b). 
— Non-negativity  of  7j  : 

From  the  definition  of  I\  by  applying  lemma  3.1  we  have. 


and 


■  Sc}\{ai?rir-  L{!r[lt]frdr- 


Also. 


1  d2y  1  d2v 


v 


r  dz 2  r  v  dz 2 


so 


2  f  1  d2v  ,  f  ,  d  1  dv  .  2  ,  j  -9  /  /  , 1  dv .  : 

4/c  /  vdr>-  (--{-—))  rdr  -  4k  /  (  — (-  —  ))rrfr  . 

Jo  r  dr^  7r«  T  dz  Jc «  r  <?r 


where  we  have  applied  (3.4c)  from  lemma  3.3. 

Therefore. 

•>  _ o  7  ,  d  1  du  2 

7 1  >  (1-4/cV,  )/  (^-(-~-))  rdr  * 

Jo  dr  r  dr 


(1 


4k2Mi  2  -  4/c4Mj  2)  j 


9  I  1  dvU2  A 
Or  r  Or  '*■ 


Since  we  want 


7)  >  0  .  by  setting  1  -  4k2^,  '  -  4k4M!  2  =  0 
=  .895  and  with  this  choice  of  k  we  have  7j  >  0  . 


we  may  pick 


K  ' 


To  bound  the  integral  term  72  ,  following  5  ,  we  may  split  72  as 


£2  ~  J  J - vJ2vdrdz  , 

£z-  f  [  -  ~^-wJ2vdrdz 

J  JRs  r‘  dr 

By  taking 

f  —  dv  dz  in  (2.4a)  we  have  that 

£  ; 

(- 

where  max  l(J, 

j  u  =  -Q  .  We  can  bound  £2  by 

C,  <  —  /  1  -  2-( J2vf)drdz  , 

*  J  JR*  r  4£r 

where  we  have 

used  the  inequality 

2  6 2 

a  6  <  sa2  -  — 

4s 

(- 

Moreover. 

[  [  1,0  f  f  \  d  \  dd  d2v  2  ,  , 

J  Jr,  t  J  J Hs  t  dr  r  dr  dz- 

f  f  t  d  1  dc-  •>  <9  1  dc  .2 1  , 

2/  jj{ar{r^]]  1 rdri:  ■ 

(- 

By  taking  £  - 

1  and  applying  (3.4c)  from  lemma  3.3.  we  have  the  following  bound  for 

r;  ,  [  I  (^('^))W-  - 

7T  y  dr  r  dr 

(- 

To  bound  £3  we  have  that  an  application  of  Schwarz's  inequality  give  us 


APPENDIX  D 


In  this  appendix  we  derive  the  expression  (4.1)  and  the  bound  (5.4). 
First  we  define  the  functional  for  s  0  by. 


S  j(tM) 


JJJ 


d  (  1  dv 
dr  r  dz 


d  1  dr)  d  1  dv  d  1  dr) 

Sb  id-B-d-rFdad-rT,' 


d  ,  1  dv  v  d  ,  1  dr)  t  d  1  dv  d  1  dr)  s  1  dv  dr)  ^ 


rdrd: 


(B -1) 


dr  '  r  dr  '  dr'  r  dr  1  dz'  r  dr  '  dz'  r  dr^  rA  dz  dz 
where  v  and  f)  belong  to  the  class  C4(B_,)  and  satisfy  the  boundary  conditions  (3. 2d),  and 
(3.2e)  and  the  regularity  conditions  (3.2/),  (3.29).  and  (3.2/i).  We  first  prove  that  the 
functional  Ss  satisfies  the  following  relation. 


a .^  =  jj  luw:- 1 


\  2  dv  d2fj  2  dv  d2ij  d  1  -  d2r)  . 

’ - -  \dr  ■ 


dr  drdz 


dz  r  dz- 


(B.2) 


We  can  prove  (B.2)  as  follows.  Starting  with  (B.l)  an  integration  by  parts  in  r  and  z 


gives, 


n  f  [  f  1  dvtd  d  ldrj 


)) 


d  d>  ld»7 

dz  dz  r  dz 


1  dv 
r  dr 


d  d  A  dr)  d  d  ldi) 

■z-  rr  -  T  T’B 

dr  dr  r  dr  dz  dz  r  dr 


\drdz 


I  (1 

JCs  T 


1  dv  d2f)  1  dv  d 2 


dz  dz2  r  dr  drdz 
Then  another  integration  by  parts  gives, 


n  ,,  f  f  1  dv  dr) 

Vz)dr  -  J  JK,r^.Fzdrdz- 


B  ,(v\  v) 


IL 


dAd.  d  l  dr)  ,  d  id  d  1  df)  ,, 

dr  r  dr  dr  r  dr  dr  r  dz  dz  r  dr 


Of  Id,  d  1 


1  d3i) 

drdz  —  / 

La.tdAdrdz.  ( 

rdrX  dr  r  dz 

r  dz3. 

J  J 

Rs  r3  dz  dz 

Using  the  definition  of  J‘  in  (3.2b)  the  first  two  terms  on  the  right  hand  side  of  (5.3) 
can  be  rewritten  as 

■  f  d  3  (  3  (l3^\\\  3  3  I  3  ll  3* \\\)  -  :  3  ll  3  t2-\  (Ua\ 

V  Xdr^rdr^dr^r  dr^  ~  dr^rdz^dzK  dr^  i  V  dr'  r  drJ  ^ 

the  next  two  terms  can  be  rewritten  as 


(£.(£  ®.(r£.(i*2)) 

^  1  d '  .1 

\  =:  VJ  \ 

\\d-  j2  _  1^] 

l  IRiM 

\dz[rdr[  dr  r  dz 

1  rdz*V  J 

I  W\ 

1  r  ds2  r3  dr2  J 

1 

the  last  term  in  (5.3)  becomes 

1  dtv  df) 
^  Yz  dz 


d  1  -dr)  1  - d2f) 

dz  r3'~  dz  r3  U  dz2 


the  two  terms  in  the  next  to  the  last  integral  of  (5.3)  can  be  expressed  as 


i-Zi  =  ±(‘ ctii -  ifsti 

~  dz 3  dzr  dz 2  r  dz  dz2 


r-d.  d  ,1  df)  v  dr)  d  d2rj  1  dv  d2r) 


r  dr  dr  r  dz  r3  dz  dr  r  drdz  r  dr  drdz 
Using  the  expressions  (5. -4a)  to  (5. 4e)  in  (5.3)  and  integrating  by  parts  we  obtain 
(5.2).  Since  £(s)  =  27rSJ(ti-,  v)  .  where  tv  is  a  solution  to  (3.2).  then  (4.1)  follows  replacing 
0  and  r)  by  v  in  (5.2)  . 
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In  the  case  that  rp  =  xp  ,  f)  =  rj  and  s  =  0  .  where  v  and  rj  are  solutions  to  (3.2)  and 
(5.1)  respectively,  we  can  prove  the  following  identity, -which  is  useful  to  furnish  a  bound 


to  £( 0), 


B  =  Bo{t},T])  ~  3q(v  -  T),W  -  T}) 


(B.5) 


which  follows  provided  that  Bo{ii',v)  —  Bo(r),r])  .  To  prove  this  last  identity  we  use  that 
Jat)  —  0  in  Rq  and  dxbjdz  —  dr)Jdz  .  dx p  j  dr  =  drj/dr  and  xb  =  77  at  z  =  0  ,  therefore  by 
using  (£.2)  we  have  that 

0,1  \  [  [  1  t4  j  j  [  (2\dT)  d2v  dr]  d'r]  ]  d  1  d2r) 

o(*'")  =  i  Jr,'"  "  !cS'  ir,37--3-r37a-J-F,[-rna^] 


Bq(t],  T]) 


Our  goal  is  to  give  a  bound  to  .£(0)  in  terms  of  £00(0),  to  do  that  we  proceed  as 
follows.  Since  rj,xj:  and  their  first  derivatives  have  the  same  value  at  the  inflow  boundary 
then,  by  using  (B. 2),  we  have 

B  o(v  -  ri,xv  -  r))  =  [  f  (- - -)J4(v  ~  v)drdz 

J  Jr „  r 


1  [  [,  <u  d  2.  1  r dxb  dJ2xb  dv  d  ,  1  2  du-  ,  t 

=  lj  jR}*-n]\-rTzJ  V~7Ai;^-TzTrJ  7-T;J  V\dri=-  (B'6’ 

where  we  have  used  that  v  satisfies  equation  (3.2 a)  and  JArj  =  0  .  Using  the  values  of  v 
and  rj  at  the  boundary  and  the  regularity  conditions  we  can  rewrite  B^{v  -  q.v  ~  77)  as 


B  o{xb  -  ti,v  -  tj)  =  ~{C  1  -  C2  ~  C3) 


(B.7) 


C'=J  ■ 

C,-// 

J  JRn  r-  dr  dz 


mIUrh&l*-"i£ 


2  (t/.’  -  r?)  d2c  y 2  1  j  j 

J  V - - — -x—r~J  u  drdz  . 

r*  oroz  J 


(v  -  rj)  r  dv  dJ2u'  2  r,  ' 

- - -  -- - - - ~—J~v',drdz 

r2  L  dz  dr  r3  dz 


(V  -  »?)  ,2  1  j  j 

- o - Vi  drdz  . 

r*  drdz  J 


We  define  5i,  S2  by 


=  IIJ. 


1  d(v  -  n)  dr\ 

^  d^  d7 


J2vdrdz 


s,m~IL?T,w 


dr?  , 

n)  —  J~vdrdz  . 

dr 


Then,  by  a  simple  algebraic  manipulation,  the  following  relation  is  satisfied 


C2-C3  =  S,  -  5: 


On  applying  Schwarz's  inequality  twice  and  (.4.5)  have  that. 


Si  < 

(/  j„J-rTr{v-,'))4rdrdZY''(j  \2S0W.v) 


a 


fl 


-A  A  V*  .m  J“  J*  J" 


S  2  < 


(/  jRii{7Tz{v-,,)]',drdz)'  \J /s„<rL5;)V"--)''4(2a»(^))‘  J' 

To  furnish  the  bounds  for  Si  and  So  we  need  an  additional  inequality 

l  f'Tdr  -2“'7(f0  (if}2rdrY  ■ 

for  /  £  C1  I (0, 1)]  with  /( 1)  =  0  . 


(5.10) 


This  inequality  follows  in  the  same  manner  as  3.4b  in  lemma  3.2  ,  Using  inequality  (5.10) 


when  f  =  \  j-r{xi>  -  rj)  ,  *£(tf  -  q)  ,  ±  §?,  and  ±  g  ,  we  get 


S  i  -  S2  <  2/ij  1(80(^  -  -  T7))1  2(So(r?.r?))1/2(S0(c'.c))1  2 


(5.11) 


Using  Schwarz’s  inequality  and  inequalities  (3.4a)  and  (.-4.5)  we  get  the  following 


bound  for  C 1  . 


Cl  <  ~~^Ml"1(^o(C''  —  7?,U>  —  77))1/2(8o(U-,  li>))  1,2 


(5.12) 


Therefore  using  these  bounds  for  C 1,  and  C2  —  C3  in  (5.11)  we  have  that 

(  Bo(*  -  r?.V  -  r,))1/2  <  l(2yl9jf±!  *  2Mr1So(r?.M)  =  )(So(v<t))1'2  . 

1/  ,T 


8  o(f  -  -  ri)  <  (2v  29.^1_)-  -  4M]  230(i7.r?)  |8.;,(f.  f)  . 

7T  J 


(5.13) 


(5.14) 


Finally  using  (5.3)  we  have  that 


5(0)  < 


£•^(0) 


js  (2\2Q^r‘  rr)  *  -  2ni~Eyc(0)  rr\ 


(5.15) 


where  we  have  used  (5.5)  and  that  B.:(rj.ri)  ~  5^(0)  2~  .  and  B.^(v.v)  =  5(0)  2~.  Then 
(5.15)  gi  ves  the  desired  bound  (5.4)  for  5(0)  in  terms  of  5^(0). 


••  .  ■  j  * 


^  V  V 


-  -"  .»  .*  ”  »  ">  »  jO 


WWW  VX 


APPENDIX  C 


In  this  appendix  we  give  a  brief  description  of  a  minimun  principle  associated  with 


the  Stokes  problem  (5.1).  We  consider  the  following  variational  problem  , 


I  ( <t>‘ )  =  min  /(<?) 

<t> 


where  I{6)  —  3o{d>,d>).  The  admissible  functions  o  belong  to  the  class  C4(R,-<)  and  satisfy 


boundary  conditions  (5.1b).  (5.1c)  and  the  regularity  conditions  (  5. Id),  (o.le).  (5. If). 


We  assume  that  the  minimun  in  (C.l)  is  attained  for  some  o' .  By  representing  o  as 


o  =  <?'  —  e\  we  consider  the  variation  SI  .  defined  by 


61  = 


dl(  O'  —  eX) 


(C.2) 


Substituting  o  by  o'  -  sx  in  H<t>)  we  have  that 


6,  -*/  J  ijVdr*  -* 


2  dx  d2©’  d  ,  X  d2©’  x  -X  do'  t 
r  dr  drdr  dc  r  dr2  r3  dc  J 


where  we  have  used  (B.2).  Since  we  want  that  61  =  0  then 


drdc  =  0 


since  x  •  and  vanish  at  z  =  0  .  Since  \  is  arbitrary  we  then  conclude  that  JAo'  =  0. 


Therefore  Eoo(0)  =  27r/(o')  . 


•r.  -\-\-\  -  -  -  .  •  j-  .-  »•  .  . 


References 


[l]  Ache,  G.  A.,  “Incompressible  Viscous  Steady  Flow  in  Channel  and  Pipes",  Ph.D 
dissertation,  The  University  of  Wisconsin-Nladison  (1985). 

[2;  Agmon.  S.,  “  Lectures  on  Exponential  Decay  of  Solutions  of  Second-Order  Elliptic 
Equations".  Mathematical  Notes,  29.  Princeton  University  Press  (1982). 

[3]  Amick,  C.  J.,  “Steady  Solutions  of  the  Navier-Stokes  Equations  in  Unbounded  Chan¬ 
nels  and  Pipes",  Ann.  Scu.  Norm.  Sup..  Pisa  4  (1977).  pp.  473-513. 

[4]  Amick,  C.  J  ,  “  Properties  of  Steady  Navier-Stokes  Solutions  for  Certain  Unbounded 
Channels  and  Pipes",  J.  Nonl.  Anal.  Th.  Appl.  6  (1978),  pp.  698-720. 

[5j  Horgan.  C.  O.,  “Plane  Entry  Flow  and  Energy  Estimates  for  the  Navier-Stokes  Equa¬ 
tions",  Arch.  Rational  Mech.  Anal.  68  (1978),  pp.  359  381. 

[6(  Horgan.  C.  O.  and  Knowles,  J.  K.,  “Recent  Developments  Concerning  Saint-Yenant's 
Principle",  in  Advances  in  Applied  Mechanics  (J.  \\ .  Hutchinson,  ed.),  23  (1983),  pp. 
179-269. 

[7]  Horgan,  C.  O.  and  Payne.  L.  E.,  “Decay  Estimates  for  Second-Order  Quasilinear 
Partial  Differential  Equations".  Advances  in  Applied  Mathematics  5  (1984),  pp.  309- 
332. 

|8j  Horgan.  C.  O.  and  Wheeler,  L.  T.,  “Spatial  Decay  Estimates  for  Navier-Stokes  Equa¬ 
tions  with  Applications  to  Entry  Flow",  SIAM  J.  Appl.  Math.  35  (1978),  pp.  97-116. 

9’  Isifyan.  G.  A..  “Saint-Venant’s  Principle  for  Flow  of  a  Viscous  Incompressible  Fluid 
and  its  Applications".  Trans.  Moscow  Math.  Soc.  J  (1983).  pp.  1-37. 


32 


[10  Joseph  D.  I).  and  Carrni,  S.,  “Stability  of  Poiseuilie  Flow  in  Pipes,  Annuli,  and  Chan¬ 


nels”,  Q.  Appl.  Math.  26  (1969),  pp.  575-599. 

[11  Knowles.  J.  K..  “On  Saint-Yenant's  Principle  in  Two-Dimensional  Linear  Theory  of 
Elasticity”  Arch.  Rational  Mech.  Anal.  21  (1966).  pp.  1-22. 

1 12;  Orr,  W.,  “The  Stability  or  Instability  of  Steady  Motions  of  a  Liquid,  Part  II.  A  Viscous 
Liquid”,  Proc.  Roy.  Irish  Acad.  Sect.  A  27  (1907),  pp.  69-13S. 

1 1 3  Skoinikoff.  I.  S..  “Mathematical  Theory  of  Elasticity” .  McGraw-Hill,  New  York  (1956). 


(14  Toupin.  R.  A..  “Saint-Venant's  Principle”.  Arch.  Rational  Mech.  Anal.  18  (1965). 
pp.  83-96. 


r'S 


